Field induced spin reorientation and giant spin-lattice coupling in EuFe2As2 
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We have studied a EuFe2As2 single crystal by neutron diffraction under magnetic fields up to 3.5 T and tem- 
peratures down to 2 K. A field induced spin reorientation is observed in the presence of a magnetic field along 
both the a and c axes, respectively. Above critical field, the ground state antiferromagnetic configuration of 
Eu^^ moments transforms into a ferromagnetic structure with moments along the applied field direction. The 
magnetic phase diagram for Eu magnetic sublattice in EuFe2As2 is presented. A considerable strain (~0.9%) is 
induced by the magnetic field, caused by the realignment of the twinning structure. Furthermore, the realign- 
ment of the twinning structure is found to be reversible with the rebound of magnetic field, which suggested the 
existence of magnetic shape-memory eff'ect. The Eu moment ordering exhibits close relationship with the twin- 
ning structure. We argue that the Zeeman energy in combined with magnetic anisotropy energy is responsible 
for the observed spin-lattice coupling. 



PACS numbers: 74.70.Xa, 75.25.-j, 75.30.Kz, 75.80.+q 

The recent discovery of iron pnictide superconductors has 
triggered extensive research on their physical properties and 
mechanism of high temperature superconductors 12-0]. All 
iron pnictides are found to be of layered structure in nature. 
For undoped iron pnictides, the chains of parallel Fe spins 
within the FeAs layers couple antiferromagnetically in the ab 
plane of the orthorhombic lattice with an antiparallel arrange- 
ment along the c axis |4-6]. This antiferromagnetic (AFM) 
order in the parent compounds is likely due to a spin-density- 
wave (SDW) instability caused by Fermi surface nesting |7]. 
Similar to the high Tc cuprate superconductors, the undoped 
iron pnictides are not superconducting under ambient pressure 
and show an antiferromagnetic SDW order. Upon carrier dop- 
ing, the magnetic order is suppressed and superconductivity 
emerges concomitantly isl^]. 

EuFe2As2 is a peculiar member of the iron arsenide 
AFe2As2 family since the A site is occupied by Eu^^, which 
is an 5"- state (orbital angular momentum L = 0) rare-earth 
ion possessing a Af structure with the total electron spin S = 
7/2. The theoretical effective magnetic moment of Eu^^ ion is 
1.94 j^B- Our previous neutron diffraction work on EuFe2As2 
single crystals revealed that the Fe and Eu spins form long 
range AFM order below 190 and 19 K, respectively fid]. Fur- 
thermore, early studies have shown that an external magnetic 
field may induce ferromagnetic (FM) order in EuFe2 As2 
which suggested a weak AFM coupling between Eu spins. 
The superconductivity can be achieved by applying high pres- 
sure or doping on either Eu or As site in EuFe2As2 il2l - 
Interestingly, some works show that the ordering of the 



Eu^^ moments persisted in the superconducting phase, such 
as in the pressure-induced EuFe2As2 superconductor [12] and 
EuFe2(Aso.7Po.3)2 superconductor [13]. However, the role of 
Eu magnetism in superconductivity is not yet clear and the Eu 
ordering state in different systems still need to be clarified. 

Here we report a single crystal neutron diffraction measure- 
ment on EuFe2As2 under a magnetic field up to 3.5 T. The 
spin reorientation of Eu moments is observed upon an applied 
magnetic field parallel to both a and c axes of the orthorhom- 
bic structure, while the Fe SDW order persists at high mag- 
netic fields. Interestingly, the application of a magnetic field 
changes the twinning population in EuFe2 As2 and the redistri- 
bution of the domain population is found to be associated with 
the evolution of the magnetic order of Eu moments, which in- 
dicates the existence of a giant spin-lattice coupling effect. A 
single crystal of EuFe2 As2 was grown by the Sn-flux method 
[10]. It was in shape of a platelet with approximate dimen- 
sions of 5 X 5 X 1 mm^. Single crystal neutron scattering mea- 
surements were performed on the thermal neutron two-axis 
diffractometer D23 at the Institut Laue Langevin (Grenoble, 
France). To investigate the evolution of the magnetic order of 
EuFe2As2 under magnetic field, the crystal was aligned in the 
a-c scattering plane and a horizontal magnetic field up to 3.5 
T was applied by using a cryomagnet. The measurement was 
performed with an incident neutron wavelength of 1.28 A. 

In zero field, the crystal structure of EuFe2As2 can be well 
described within the orthorhombic symmetry at 2 K and the 
magnetic reflections originated from both Eu and Fe magnetic 
sublattices were observed. The long range ordering of Eu^^ 
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FIG. 1: (Color online) (a) Evolution of (203)m magnetic reflections 
with the change of magnetic field at 2 K.(b) Magnetic field depen- 
dence of (203 )m magnetic reflection in two dimensional plot, (c) 
Evolution of (200) o reflection with the change of magnetic field at 
2 K. (d) Magnetic field dependence of (200)<9 reflections in two di- 
mensional plot. 



[with k = (0,0,1)] and Fe^+ moments [with k = (1,0,1)] forms 
a AFM structure as confirmed in our previous neutron diflTrac- 
tion study [10]. Several typical AFM reflections of the Eu 
ordering, such as (201)m, (203)m and (005)m, were selected 
to examine the magnetic structure evolution under the mag- 
netic field {H) along both c and a axes, i.e. [001] and [100] 
directions, respectively. As an example, the field dependence 
of the (203 )m reflection is plotted in Fig. 1(a) and Fig. 1(b) in 
two dimensions. Application of an field along c axis strongly 
weakens the (203 )m reflection and suppresses it totally at the 
critical field (if^;?.^) of 1.8 T at 2 K. On the other hand, the 
increase of intensity at some nuclear reflection positions, e.g. 
(200)(9 as shown in Fig. 1(c) and (d), indicates that the Eu spin 
gradually reorientate from the ab plane to the [001] direction. 
The field-induced magnetic phase transition takes place from 
antiferromagnetic, via a canted configuration, to the ferromag- 
netic structure. Note that the small nuclear structure factor of 
the (200) o reflection makes it possible for us to clearly ob- 
serve the contribution from magnetic scattering. 

Fig. 2. (a) and (b) display the field dependence of inte- 
grated intensities of selected magnetic scattering reflections 
in EuFe2As2 with the field parallel to the c axis at diflTerent 
temperatures. The critical field at which the Eu magnetic or- 
der changes from AFM to EM is decreasing with increasing 
temperature. When the field is parallel to the a axis, it is 
observed that the intensity of the (201)m reflection increases 
slightly and then decrease sharply with increasing strength of 
the applied field at 2 K [see Fig. 2(c)]. It is known that diflfrac- 
tion signals from both (201)m and (021)m reflections will be 
detected due to the twinning configuration [10]. The initial 
increase of the (201)m intensity in Fig. 2(c) is caused by the 
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FIG. 2: (Color online) (a)(b) Integrated intensity of (203)m and 
(200) <9 reflections as a function of magnetic field with the field ap- 
plied parallel to c axis. (c)(d) Magnetic field dependence of inte- 
grated intensity of (201)m and (006) o reflections at difi'erent tem- 
peratures with the magnetic field applied parallel to a axis. (e)(f) 
Temperature dependence of integrated intensity of (200) o and (006) o 
reflections. 

increasing contribution from the (021)m reflection, which is 
related to the rotation of the (0^0) domains. Actually, the re- 
distribution of domains under an applied field is observed and 
it will be discussed in the following text. In Fig. 2(d), the 
increase in the intensities of the (00/) reflections with / even 
is observed instead of the (/zOO) reflections. All these results 
suggest that the Eu spins will orient to the direction of the 
applied magnetic field once the field strength is greater than 
the critical field. The EM arrangement of Eu spins is further 
confirmed by the refinement of the reflections collected un- 
der the field of 3 T at 2 K. The amplitude of Eu moment is 
estimated to be 6.9(4) yug and 6.7(4) yug for the two EM struc- 
tures with Eu moments along the c and a axis, respectively. 
The temperature dependence of EM reflections was measured 
afterwards to determine the ordering temperature of the Eu 
moments. Fig. 2(e) and (f) show the temperature evolution of 
integrated intensities of the (200)o and (006)^ Bragg reflec- 
tions for H parallel to c and a axes, respectively. A power-law 
was adopted to describe this second order transition by fitting 
the temperature variation in the ordering parameters. As in- 
dicated by the dash line, the power-law fit on the data at 3 T 
in Fig. 2.(e) yields a magnetic ordering transition tempera- 
ture Tc = 33.5(4) K. The intensity tail above Tc is caused by 
the contribution from the field induced moment. According to 
the neutron measurement results, the magnetic phase transi- 
tion temperatures are determined and plotted in Fig. 3(a) and 
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FIG. 3: (Color online) (a) Magnetic phase diagram for EuFe2As2 with applied field parallel to the crystallographic c axis. The inset shows the 
schematic view of the magnetic structure. The AFM magnetic structure at zero field is denoted as the open arrows, while the field induced FM 
structure is denoted as the solid one. (b) Magnetic phase diagram for EuFe2As2 with applied field parallel to the crystallographic a axis. The 
lines are guide to the eye. See text for more details. 



(b). It can be seen that the application of a magnetic field not 
only changes the ordering configuration but also enhances the 
ordering temperature of the Eu moment. 

Besides of Eu order, we also examined the Ee order un- 
der the applied field along both c and a axes. In contrast 
to the reorientation of the Eu spins, the AEM SDW order 
of Ee is found to be robust and it persists till fields up to 
3 T. The two dimensional plot of the (101)m reflection un- 
der a 3 T field at 2 K is shown in Eig. 4(a). Surprisingly, 
it is found that the integrated intensity of the (101)m mag- 
netic reflection decreases slightly and then sharply increases 
with increasing magnetic field, as shown in Eig. 4(b). As we 
known that a twinning structure may exist in orthorhombic 
EuEe2 As2 phase due to the interchange of the orthorhombic a 
and h axes. Additionally, the intensity of the (101)m reflec- 
tion only originates from the magnetic scattering of the (/zOO) 
twin. In order to clarify the evolution of twinning structure 
under applied field, the (400) o reflection was examined from 
zero field up to 3.5 T since both (040)o and (400)^ reflec- 
tions measured in present magnetic structural configuration 
are pure nuclear reflections. As shown in Eig. 4(c), the twin- 
ning structure forms at zero field and results in the splitting 
between the (400)o and (040)^ reflections. The domain pop- 
ulation is estimated to be 1:3 for the (/zOO) and (0^0) twins 
at zero field. It is interesting that only the (400)^ reflection 
is observed when the magnetic field is increased up to 3 T 
[Eig. 4(d)]. The single reflection from (400)o is further con- 
firmed by carefully checking the vicinity in reciprocal space, 
which indicates that the single crystal is detwinned. The de- 
tailed evolution of the domain population with increasing field 
is illustrated by the intensity change of the (400) o and (040) o 
reflections [Eig. 4(e)]. With increasing field, the intensity of 
the (400) o decreases slowly and reaches the lowest at the crit- 
ical field H^li^ and suddenly increases at the higher critical 



field H^^.^. Whereas the (040)o reflection vanishes at H^^.^ 
in the other hand. Therefore, H^^.^ can be considered as the 
critical field at which the (0^0) domains overcome the domain 
wall energy and realign to (/zOO) domains. The realignment 
of the (0/:0) domains will firstly introduce internal stress and 
this stress will act on the (/zOO) domains and tilt them slightly 
away from the balance position, which leads to the decrease of 
the (400)o reflection. However, (/zOO) domains will rebound 
to the balance position when the (0^0) domains realign to the 
(/zOO) domains. Therefore, H^l^ corresponds to the critical 
field when the (/zOO) domains tilted and formed the largest tilt- 
ing angle away from the balance position. The Q-scans of the 
(400)o reflection under three typical fields are plotted in Eig. 
4(f). Based on the peak position of the (400)o and (040)o re- 
flections, the strain caused by the twin boundary motion was 
estimated to be around 0.9%. In contrast with some magnetic 
shape-memory alloy, where the field induced strains are asso- 
ciated with the martensite/austenite phase transition |[l6[[l7 1. 
the field induced strains in the present EuEe2As2 single crys- 
tal is caused by the domain realignment purely. Furthermore, 
it is observed that the realignment of domains is a reversible 
process with a rebound field, as shown in Eig. 4.(c). 

The critical fields H^l^ and if ^^.^ are also determined at dif- 
ferent temperatures and plotted in Eig. 3 as the red circles 
and dots. It is interesting that the critical field of the domain 
redistribution is closely correlated with the field that induces 
the AEM-EM transition of the Eu magnetic sublattice. The 
energy diflTerence between two domains seems strongly re- 
lated with the total energy of the Eu magnetism. Usually, the 
application of mechanical compressive stress a on twinned 
single crystals can induce a motion of the twin boundary, so 
called 'detwinning' process. For a twinned single crystal with 
EM structure, the application of a magnetic field can also 
introduce stress which is generated through the Zeeman en- 
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FIG. 4: (Color online) Selected Bragg reflections of EuFe2As2 at 2 
K. (a) The contour map of (101)m reflection at 3 T field, (b) Field de- 
pendence of integrated intensities of (400)^ and (101)m reflections, 
(c) and (d) The contour maps show the Q dependence of the (400)o 
and (040)o nuclear reflections at zero field and 3 T field, respectively, 
(e) Evolution of the distribution of (400)<9 and (040)<9 reflections with 
increasing and decreasing magnetic field at. (f) g- scans of (400)<9 re- 
flection under three typical magnitudes of magnetic field. 



ergy E = -jaqM-H and magnetic anisotropy energy. This de- 
twinning process can be realized if the field induced stress is 
larger than the internal stress originated from twin disloca- 
tion. As illustrated in Fig. 3(b), the Eu spin rotates toward 
the a axis with increasing field. Concomitantly, the Zeeman 
energy decreases with increasing field strength and magneti- 
zation of the EuFe2As2 crystal. Thus, by minimizing the total 
energy, greater stress is generated and results in a single do- 
main. Since the twinned structure only exist in the ah plane in 
present orthorhombic structure, the domain population does 
not change when a magnetic field is applied parallel to the 
c axis. Note that the twinned structure also plays a key role 
with respect to the magnetic and electric behaviors. For exam- 
ple, the in-plane resistivity can be aflTected remarkably by twin 
boundary scattering and the magnetization is associated with 
the response of diflTerent domains to the applied field. Since 
the twinning structure exhibits as the common feature for al- 
most all orthorhombic phases in iron pnictides It IS 
of great significant to detwin the crystals before performing 
various experiments to evaluate the real intrinsic properties. 

Within the frame of the standard model for rare earth, the 
crystalline electric field (CEF) eflfect is responsible for the 
magnetic anisotropy for most of the rare earth ions with finite 



orbital magnetic moment. However, for 5"- state rare earth ions, 
such as Gd^^ and Eu^^, the CEF eflfect is negligible because of 
the vanishing of the orbital magnetic moment and the charge 
density is no longer coupled to the spin. Nevertheless, weak 
magnetic anisotropy is observed in both Gd^^ and Eu^^ con- 
tained compounds and the magnetic anisotropy is argued to 
be driven mainly by dipole-dipole interactions 1 2Ql 422ll . The 
magnetic anisotropy energy caused by the dipole-dipole inter- 
actions in EuFe2As2 was evaluated for three diflferent ordered 
states of Eu^^ magnetic moments: ground state AFM config- 
uration, FM configurations with moments aligned along a and 
c axes, respectively. By considering the contribution of neigh- 
boring magnetic atoms within the sphere of 25 A radius, the 
dipolar energies are obtained to be -206 yueV/Eu, -72 yueV/Eu 
and 158 yueV/Eu for above mentioned three ordered configu- 
rations, respectively. It clearly suggested that the dipolar in- 
teraction favors the AFM ordered state, while the hard axis 
is predicted to be c axis. Our experimental data are in good 
agreement with the prediction from the dipole-dipole interac- 
tion, which indicates a dominant contribution of the dipole- 
dipole interactions to the magnetic anisotropy of EuFe2As2. 
Given the fact that the magnetization saturated at diflTerent 
fields for H parallel to the a and c axis due to the magnetic 
anisotropy rUD, we argue that the magnetic anisotropy en- 
ergy may also generate stress and be partly responsible for the 
movement of the twin boundaries under field. 

In summary, our single crystal neutron diflTraction experi- 
ments show a magnetic field induced magnetic phase transi- 
tion in EuFe2As2 with the Eu moments changing from AFM 
to FM arrangement. The ordering temperature of the Eu mo- 
ments increases with increasing field. Moreover, giant spin- 
lattice coupling has been observed as indicated by the redis- 
tribution of the domain population. Since the domain realign- 
ment is intimately correlated with the magnetic phase transi- 
tion, the spin-lattice coupling in EuFe2As2 can be attributed 
to the stresses generated by Zeeman energy and magnetic 
anisotropy energy. 
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